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In order to construct a general density-functional theory for nonuniform fluid mixtures, we 
propose an extension to multicomponent systems of the weighted-density approximation (WDA) of 
Curtin and Ashcroft [Phys. Rev. A 32, 2909 (1985)]. This extension corrects a deficiency in a 
similar extension proposed earlier by Denton and Ashcroft [Phys. Rev. A 42, 7312 (1990)], in that 
that functional cannot be applied to the multi-component nonuniform fluid systems with spatially 
varying composition, such as solid-fluid interfaces. As a test of the accuracy of our new functional, 
we apply it to the calculation of the freezing phase diagram of a binary hard-sphere fluid, and 
compare the results to simulation and the Denton- Ashcroft extension. 
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Density-functional theory has been applied with suc- 
cess to the calculation of the equilibrium properties of a 
variety of inhomogeneous fluid systems H. One of the 
most successful of these theories has been the Weighted- 
Density Approximation (WDA) of Curtin and Ashcroft 
||. The WDA has been applied with good results to 
problems of solid-fluid phase equilibria |^J3], solid-liquid 
interfaces anc ^ nonuniform fluids For sys- 

tems with well-defined bulk densities that are uniform 
throughout (such as in the case of freezing calculations), 
the Modified WDA (MWDA) of Denton and Ashcroft 
jjj-H may be used, yielding essentially identical results 
as WDA with much less computation. By construction, 
however, the MWDA cannot be used to study systems 
in which the bulk densities vary over length scales larger 
than a few atomic diameters (for example, interfaces or 
nucleation), so the full WDA must be utilized in such 
cases. (It is, however, often possible to exploit the system 
symmetry to allow for some simplification of the WDA 
scheme, e.g. the Planar WDA of Marr and Gast Jio[.) 

Denton and Ashcroft have proposed an extension of 
the WDA and MWDA for multi-component systems jll| . 
Using the extended MWDA they calculated portions of 
the phase diagram for the binary hard-sphere system and 
obtained excellent agreement with the simulation results 
of Kranendonk and Frenkel jl2). However, the WDA 
that they proposed assumes that the mole fractions of 
the various components of the mixture are global prop- 
erties of the system, making their approach inapplicable 
to systems, such as interfaces, in which spatially extended 
variations of composition may exist. In this work, we pro- 
pose a different extension of WDA for multi-component 
systems that correctly preserves the local nature of the 
WDA and thus can be applied to any inhomogeneous 
fluid system. In order to verify the validity of our ap- 
proximation, we also formulate the extension of MWDA 
consistent with our extension of WDA and compare our 
results for freezing of binary hard-sphere mixture with 
those of Denton and Ashcroft Jll| . 

We formulate our extension of WDA for the case of 
binary mixture - the extension to three or more com- 
ponents is straightforward. The state of a binary mix- 



ture is specified by the single-particle densities of each 
species pi(r) and P2(r), or, alternatively, by the total 
density p(r) = pi(r) + ,02(1") and the mole fraction x(r), 
which we take to be the concentration of component 
2 ( x = X2 — p%/p ). A fundamental result of the 
density-functional theory is that the Helmholtz free en- 
ergy J-[pi,p2\ is a unique functional of the densities and 
is minimized (at constant average density) by the correct 
equilibrium densities. 

The free energy can be separated into ideal and excess 
terms as follows 



F[Pl ,P2\= F\& [Pi , P2\ + Fox [Pi , P2\ 



(1) 



where the ideal part is simply the sum of the ideal free 
energies of the individual components 



^id[pi,p 2 ] = P^Y,J drpi(r)[lnA? Pi (r) - 1] 



(2) 



with A, being the thermal wavelength of component i. 
An analytic expression for the excess free energy is not 
available and, in practice, this term must be approxi- 
mated. In general, a density-functional theory is a pro- 
cedure for the specification of this excess free energy func- 
tional. One natural starting point for the development 
of such theories is the hierarchy of n-particle direct cor- 
relation functions c-™\ defined as functional derivatives 
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,j = l,2. (3) 



Most theories to date are constructed so that Eq. (||) is 
exact in the homogeneous fluid limit - at least for n = 2. 
Of such theories, the WDA and MWDA have been per- 
haps the most successful for single-component systems. 
In what follows we discuss extensions of these theories 
to multi-component systems - the original single compo- 
nent theories can be recovered in the limit that the mole 
fraction of one species becomes unity. 



1 



Denton and Ashcroft generalize the WDA to binary 
mixtures by writing the WDA excess free energy in the 
form (l| 



J r ex[px,P2]= / dr pi(r)ip(p {1) (r),x) 

+ [drp 2 (r)iP(pW(r),x), (4) 



where ip(p,x) is the excess free energy per particle of 
the corresponding uniform mixture, and x represents the 
global average concentration of the nonuniform mixture. 
The two different total weighted densities p^(r) and 
p^(r) are defined as weighted averages of the physical 
densities with respect to weight functions u>ij according 
to 

(r) = J2 I dr ' Pi - r'| ; /5« (r) , x) , 

j=i J 

i=l,2. (5) 
The weight functions must be normalized 

Jdv'w ij {\v-v'\;p) = l, i,j = 1,2, (6) 

and such that the approximate .F cx exactly satisfies 
Eq. (|J) for ?i = 2 in the uniform fluid limit 



Ctj(|ri - r 2 |;pi,p 2 ) 
= — (3 lim 



S^exlpl, pi] 



8pi{v 2 ) 8pj{r\) 



»,J'=1,2. (7) 



In order to compute the weight functions, one now has to 
solve the system of three nonlinear differential equation 

& 

The MWDA is generalized by Denton and Ashcroft in 
the same spirit with 



•fex[pi, Pi] = Ni ip(p {1 \x) + N 2 i>{pW,x) , 



,"(2) 



(8) 



where the weighted densities p^ and p^ are position- 
independent and defined as 



2 

W^J drdr ' pi ^ pj ^ 

1 3=1 



r\;p 
= 1,2 



a;) , 



(9) 



As in the single-component case, the MWDA is much 
less computationally demanding, so Denton and Ashcroft 
have used this approximation to compute freezing con- 
ditions of the binary hard-sphere mixtures with different 
diameter ratios. Their results for a > 0.85 closely follow 
those obtained in the simulations by Kranendonk and 
Frenkel gf|. 

Nevertheless, the binary mixture WDA of Denton and 
Ashcroft defined by Eqs. ([I]) and (^) cannot be used to 
study any binary system with extended spatial variations 
of average composition. The problems are best illus- 
trated by the example of a planar crystal-fluid interface. 



The excess free energy JF ox in Eq. and, therefore, the 
weight functions tUy in Eq. (||) are defined to depend 
explicitly on the global average concentration x. Obvi- 
ously this quantity cannot be uniquely defined for the 
interfacial system, where it has to be equal to the crystal 
coexistence values x c on the one side of the interface and 
to the fluid value x on the other side. Matematically, 
the WDA equations for Wij (k; p, x) derived by Denton 
and Ashcroft (ll]] contain terms proportional to 8k,o, so 
that the computed weight functions are discontinuous in 
Fourier space at k — 0. This means that one cannot re- 
cover the bulk crystal or fluid properties in the regions 
remote from the interface (in the limit z ±oo, where 
z is the axes perpendicular to the interfacial plane). 

In order to circumvent these difficulties, we pro- 
pose an extension of the WDA to binary mixtures in 
which weighted concentrations x^(r) and x^(r) are 
introduced in addition to the weighted densities. The 
weighted concentrations replace the average concentra- 
tion x in Eq. (|) and are defined according to 

*®( r ) = lM7^ /drV 2 (r') 

i = l,2. (10) 

It is more convenient to write this extension scheme 
in terms of the weighted densities of individual species 
p*j\r) defined as 



p(i)(r) 

xw i2 (\r-r'\;p^(r),x^(r)) 



pf(r) ee J dr' Pj (r') Wij (\r-r'\;pr(r),p^(r)) , 

i,j = 1,2, (11) 
which are input into the excess free energy expression 



'I. *WfV\ 



+ /drp 2 (r)V.(M 2) (r),4 2) (r)) 



(12) 



where ip(pi ,p 2 ) is the same excess free energy per particle 
as in Eq. (^), but expressed in variables p\ and p2- 

Substitution of .F ex from Eq. Q into Eq. (0) leads 
explicitly to the following set of equations for the three 
weight functions in Fourier space 
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dp 
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(13) 
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dpi dp 2 
dip 



piwn(k) + p 2 w 22 (k) w\ 2 (k) 



dpi 
dip 
dp 2 



dwu(k) dwi 2 {k) 
piw 12 (k) — 5 h p 2 w 22 ik) 



dp 2 

dwi 2 (k) 
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dp 2 
dw 22 (k) 
dpi 



(14) 
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FIG. 1. WDA weight functions Wij(k) for a binary 
hard-sphere mixture. The system parameters are: a = 0.8 , 
x = 0.7, and p — 1.1 
curves correspond to Wu . 



. The dotted, solid, and dashed 
W12 , and W22 , respectively. 
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(15) 



The normalization of the weight functions ensures that 
at k = Eqs. (|l3|)-([l5|) correctly reduce to the compress- 
ibility rules for a binary mixture [pL4j] 



(3 1 c ij (k = 0;pi,p 2 ) 



dip 
dpi 



dtp 



- P 



d 2 ip 
dpi dpj 

*,i = i,2 



(16) 



Unlike the original generalization of Denton and Ashcroft 
jyj , the above equations do not contain terms propor- 
tional to Sk,o, and are thus continuous functions of k. 

In order to calculate the weight functions from 
Eqs. (|l3|)-(|l5|), we need to specify Cjj(fc), the Fourier 
transforms of the uniform direct correlation functions, 
and ip(pi,P2)i the Helmholtz free energy per particle. 
As an example, we have calculated the weight functions 
for the binary hard-sphere mixture with diameter ratio 
a = CTi/o"2 = 0.8 and use as input the analytical expres- 
sions for a uniform binary mixture within the Percus- 
Yevick approximation |ll||l4j]. 

Eqs. (|l3|)-([l5|) is a system of nonlinear differential equa- 
tions. We use an iterative scheme, which has satisfactory 
convergence for all values of k and density. We start by 
assuming that w±2 and all derivatives of itty are zero and 
solve the quadratic equations for w±i and W22- Then we 
estimate the derivatives of Wu and W22 using finite differ- 
ences and solve equation for W12, which is linear in W\2- 
Finally, we estimate the derivative of w\2- After 5-10 it- 
erations the solution converges to within five significant 
digits. Fig. [l] shows loy (k) computed for the binary hard- 
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FIG. 2. Comparison of binary hard-sphere crystal-melt 
phase diagrams computed with the original MWDA by Den- 
ton and Ashcroft |llj] (dashed lines) and our new extension 
of the MWDA (solid lines) . Also shown are the rescaled sim- 
ulation data from Ref. (see text for explanation). 

sphere fluid with diameter ratio a = 0.8 concentration 
x = 0.7, and total density p = 1.1 o^T 3 . All three func- 
tions smoothly approach unity when k — > 0, thus satis- 
fying the normalization conditions and allowing applica- 
tion of the WDA to binary systems with spatially varying 
characteristics. 

In order to compare our extension of the WDA to mix- 
tures with that of Denton and Ashcroft [|llj, we repeat 
the freezing calculations for binary hard-sphere mixture 
using MWDA consistent with the extension of the WDA 
presented earlier. The excess free-energy functional is 
written as 



N 1 iP(p[ 1 \p ( i ) ) + N 2 ^(p 1 2 \p2 2) ), 



where the weighted densities p[^ are defined as 



(17) 



pf = i-|drdr'p 4 (r)p,(r'K 3 (|r ~r'|;p« , 

i,j= 1,2, (18) 

This extension of MWDA yields the following equations 
for the weight function Wij (k) in Fourier representation 
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dpi 



dip 
d P] 



Wij(k) + 5 k ,op 



d 2 ip 



dpidpj 
*,i = l,2 



(19) 



where p = p± + p2 is the total density. Agreement with 
the compressibility rules given by Eq. [l^ is obvious. 

Comparison between the freezing parameters obtained 
within our extension of the MWDA and that of Denton 
and Ashcroft is shown in Figs. |[^. Fig. ^ presents the 
pressure-concentration phase diagrams together with the 
simulation results of Kranendonk and Frenkel [Q , which 
are scaled down by a factor of 0.871 in order to match 
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FIG. 3. (a) Packing fraction in the crystal phase at co- 
existence. Dashed lines represent the results of the original 
MWDA and solid lines represent the results from our exten- 
sion of the MWDA. (b) Packing fraction in the fluid phase at 
coexistence. 
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FIG. 4. Lindemann parameters for the two types of spheres 
calculated according to the original MWDA (dashed lines) 
and our extension (solid lines). 

putational Science (KITCS) and the Kansas Center for 
Advanced Scientific Computing (KCASC) for support 
and computing resources. 



the transition parameters in the single-component limit 
(x = and 1). The scaling compensates for the discrep- 
ancies between theory and simulation in the freezing pa- 
rameters of the one-component system and helps to em- 
phasize the general features of the binary system phase 
diagram. Overall, the two extensions of the MWDA pro- 
duce very similar phase diagrams and reproduce fairly 
well qualitative behavior of the simulation results. 

Fig. H compares the crystal and fluid packing fractions 
at coexistence. The packing fractions in the crystal phase 
differ for the two extensions much more than the fluid 
p ackingf r act ions . 

Fig. shows the Lindemann parameters for the two 
types of spheres. In both cases, the smaller spheres have 
larger Lindemann parameters for all values of concentra- 
tion. The Lindemann parameters for the two extensions 
of MWDA coincide for larger spheres at x = 1, and for 
smaller spheres at x = 0, where the single-component 
limit is recovered for the remaining type of spheres. 

In conclusion, we have constructed an extension of 
the Weighted-Density Approximation (WDA) of Curtin 
and Ashcroft Q to multi-component systems. This new 
density-functional theory is more general than a previous 
extension by Denton and Ashcroft 1 1 in that it preserves 
the local character of the weighted densities and thus al- 
lows its application to the systems with extended spatial 
variations in average composition, such as solid-liquid in- 
terfaces. 
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